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We consider the minimal model in which dark matter is stabilized by a non-Abelian discrete 
symmetry. The symmetry group is taken to be Dz — S3, which is the smallest non-Abelian finite 
group. The minimal model contains (nontrivial) singlet and doublet scalar representations of D3 
which couple to the Standard Model fields via the Higgs portal. This construction predicts two 
species of dark matter over much of the parameter space. Nontrivial interactions under D3 lead to 
a novel thermal history of dark matter, while the multi-component nature of dark matter can be 
tested by future direct detection experiments. 
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Understanding the nature of the cosmological dark 
matter (DM) that constitutes one quarter of the energy 
density of the universe is a central goal of particle physics 
today [l[ . While there is little room left to doubt the ex- 
istence of DM, its microscopic properties are virtually 
unknown. One of the few properties in which we can 
be confident is that DM should be stable on time scales 
greater than the age of the universe, suggesting the ex- 
istence of a new "dark" symmetry. But precisely what 
symmetry stabilizes DM is a mystery. 

Many models employ a discrete Zi symmetry to sta- 
bilize DM. This Z2 symmetry is often motivated by the 
need to suppress dangerous operators in new physics sce- 
narios that solve the hierarchy problem. However, given 
that we have no experimental indication of what new 
physics, if any, addresses the naturalness issues in the 
SM, one may take a more general perspective regard- 
ing DM and the symmetries responsible for its stability. 
More pragmatically, the exploration of alternative stabi- 
lizing symmetries is warranted by the prospect of novel 
phenomena associated with DM, as such symmetries may 
predict new states and interactions. 

Indeed, there are many possibilities other than a im- 
parity that can protect DM against decay. In partic- 
ular, besides the Abelian cyclic symmetry Zn [2j, DM 
may well be stabilized by a non-Abelian discrete symme- 
try. Non-Abelian finite groups have received some lim- 
ited attention within the context of DM. Motivated by 
improved gauge coupling unification, Ref. [3( considered 
an additional Higgs doublet in a non-Abelian discrete 
multiplet serving as DM. Non-Abelian discrete symme- 
tries also lead to distinct decay patterns in decaying dark 
matter scenarios 0, [B| . Continuous non-Abelian symme- 
tries originating from broken or confined gauge theories 
can also ensure DM stability p. Models of DM stabi- 
lized by Abelian discrete symmetries that descend from 
higher non-Abelian ones have been motivated by astro- 
physical anomalies Q , discrete gauge symmetries Q , and 
neutrino physics Indeed, non-Abelian flavor symme- 
tries are widely used to explain the neutrino oscillation 
data [lot (for recent reviews see Refs. 11 , 12 1). Such non- 
Abelian discrete symmetries can come from the breaking 



of continuous flavor symmetries [13| or from orbifold com- 
pactification of extra-dimensions [14j . 

In this paper we construct the minimal model of DM 
in which stability is a consequence of a non-Abelian dis- 
crete symmetry. The model follows in spirit the canonical 
scalar Standard Model (SM) singlet models of Refs 
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17[ . The symmetry group we consider is D3 , the dihedral 
group of order 6, which is the smallest non-Abelian finite 
group. It is isomorphic to S3, the permutation symmetry 
of three objects. We determine the minimal field content 
of the model and couple it to the SM via renormaliz- 
able Higgs portal interactions. Over much of the param- 
eter space, the D 3 symmetry predicts two species of DM, 
which can contribute to the cosmological abundance. We 
investigate in detail the cosmology of the model, finding 
that nontrivial interactions under D3 in the dark sector 
lead to a novel thermal history. We also analyze the con- 
straints and prospects at direct detection experiments. 
In particular we find that the two-component nature of 
DM can be discovered in future ton-scale experiments. 

D3 stabilization of dark matter. We start by 
constructing the minimal model of DM stabilized by a 
non-Abelian discrete symmetry. The model is based on 
the group D3, which is the smallest non-Abelian discrete 
symmetry group and describes the symmetry properties 
of the equilateral triangle. The group contains two gen- 
erators, A and B, obeying the properties 



A 3 = 1, B 2 = 1, ABA = B, 



(1) 



and the six group elements are constructed through com- 
binations of the two generators: 1, A, A 2 , B, BA, BA 2 . 
The group D3 has two singlet and one doublet represen- 
tation, denoted, 1 1: 1 2 , and 2 respectively. The singlet 
\ 1 transforms trivially under D3, A = B = 1, while the 
generators of the (nontrivial) singlet representation 1 2 
are A = 1 and B = — 1. The generators for the doublet 
representation 2 are 



A = 



£ 2-Ki/3 






-27TI/3 



B = 



1 

1 



(2) 
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The multiplication rules for the singlet representations 
are given by 

Ii ® Ii = Id I2 ® I 2 = Ii> and Ii <8) 1 2 = I 2 , (3) 

whereas the product of two doublets, 2 ® 2 decomposes 
into two singlets and one doublet as 



(0162 + ^2^2) ~ Ii, (oi&2 - a 2 & 2 ) ~ I 2 



(4) 



where (ax, a 2 ) T , (61, &2) T ~ 2; for further formulae of di- 



hedral groups, see, e.g., [12 
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The minimal model of DM having nontrivial properties 
under D3 contains two scalar fields, 



V, 



X 



- I x 

X* 



(5) 



which transform as a singlet 1 2 and doublet 2, respec- 
tively. The scalar 77 is real while x is complex, which 
amounts to a total of three new degrees of freedom. 
Along with the SM Higgs field H , the scalar sector of the 
Lagrangian contains kinetic terms and a potential. The 
general renormalizable scalar potential invariant under 
D3 and the SM gauge symmetries may be constructed 
using Eqs. (|3|4[) and reads 



T/ 2zjtiT 1 1 22, 2* Ml ^ 3 *3\ 

V = m x WH + -m 2 rj + m 3 x X + -^{X + X ) 



+ A 1 (i/ttf)2 + ^ r/ 4 + A3(x * x) 2 



ai {rfH)rf + 2a 2 (H^H)( X *x) + <W(x*x) 



(6) 



where all parameters in the Lagrangian are real. We note 
that that the theory is invariant under P and C — a fact 
that will simplify our considerations regarding the relic 
abundance. 

A model with only one scalar field 77 (x) can provide a 
viable theory of DM, but such a theory is equivalent to a 
theory based on an Abelian discrete Z 2 {Z3) symmetry. 
The minimal theory based on the group D3 contains both 
scalars, 77 and Xj and the nontrivial interaction predicted 
by this symmetry is the last term in the potential ([B]), 
with coupling 04. 

We will be interested in the case in which the elec- 
troweak symmetry is spontaneously broken while the D3 
discrete symmetry is unbroken, 



s/2\v 



(V) = 0, (x) = 0, (7) 



where v 2 = — mf/Ai ~ (246 GeV) 2 is the (squared) Higgs 
vacuum expectation value. We demand that the poten- 
tial is bounded from below, and that the electroweak 



vacuum is the global minimum of the potential. These 
requirements constrain the values of the coupling con- 
stants in the potential. Qualitatively, the quartic cou- 
plings Ai,a,3 must be positive, ai j2i 3 must be greater than 
some minimum (negative) value, and the magnitude of 
q;4 must be smaller than some maximum value. Further- 
more, the mass parameters m\ 3 should not be too large 
and the magnitude of the cubic coupling [i\ is bounded 
from above. Throughout our analysis, we have verified 
numerically that our parameter choices lead to a stable 
potential and electroweak vacuum. One may also investi- 
gate the vacuum structure of the theory at higher energy 
scales using a renormalization group analysis (as, e.g., in 
Ref . [20( ) , but here we will be content to consider Eq. © 
a low energy effective theory of DM defined around the 
weak scale. 

Expanding around the vacuum ([7]) with v — > v + h, we 
derive the following potential: 

t r 1 2 1,2 1 „2 2 1 2 * 

v = 2 m h h + 2 m >j 7 7 + m xX X 

+ X lV h 3 + axvhrf + 2a 2 vh(x*x) + ^(x 3 + X* 3 ) 

+ ^ + ^ + Hx*x) 2 

+ y^V + a 2 h 2 ( X *x) + a 3 il 2 (x*x) 
+ -gj->7(X -X )> 
where the physical masses are 



(8) 



m l 


= 2Aiu 2 , 










m 2 


= m\ + a\v 2 , 




m\ 


= m\ + a 2 v 2 . 


(9) 



The nontrivial interaction predicted by the D 3 sym- 
metry, with coupling 014, allows the scalar rj to decay via 
77 — > 3x,3x* if m v > 3m x . In this case, x is the only 
stable DM candidate. If, however, m v < 3m x these de- 
cays are kinematically forbidden, and both 77 and x are 
stable. Therefore, the D3 model predicts two species of 
DM for rrirj < 3m x . 

Cosmology. We now examine the thermal history of 
the 77 and x particles in the minimal D 3 model. The cos- 
mic relic abundances of these particles are governed by 
Boltzmann equations that account for the expansion of 
the universe as well as for reactions that change the total 
number of particles for the species of interest. If these re- 
actions freeze out when the temperature is not too much 
smaller than the mass of the particle, there will be as- 
sociated a significant energy density that remains today 
that may account for the observed DM in the universe. 

From the potential in Eq. ((8]), we can determine the 
reactions that change the total number of particles, which 
may be classified as follows: 
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(a) Annihilation into SM: 

TjTj -> ti, hh, ZZ, WW, 66 ... , 
XX* -> hh, ZZ, WW, bb... . 

(b) Semi- annihilation: 

XX -> h X*, Xh -> X*X* ■ 

(c) DM conversion: 

VX^X*X*, VX* -> XX, xx->to* 
VV^XX*, XX* -> Wi- 



ld) Late decay: 



V -> 3x, 3x* 



The reactions listed above constitute the relevant pro- 
cesses that change the number of 77 and x particles; analo- 
gous reactions hold for x* ■ (Annihilation into three body 
final states can be important right below the WW thresh- 
old [UJ.) We see that the D3 symmetry allows for the 
non-canonical processes of semi-annihilation (b) and DM 
conversion (c), studied recently in Refs. (gI. |22|. as well 
as a late decay scenario (d) [23]. CP invariance of the 



potential ([6]) implies 



so that we only need to 



solve the Boltzmann equation for n x ; the total relic den- 
sity is then = 2n x . Furthermore, processes (6), 
(c), and (d) couple the Boltzmann equations for 77 and 
X- 

We have performed a general analysis of the Bol tzmann 
equations following the usual treatments of Refs. 
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with modifications for semi-annihilation and DM conver- 
sion; a summary is presented in Appendix A. In the re- 
sults presented below, we obtain the relic densities of 77 
and x by numerically integrating the coupled set of equa- 
tions. We now present a detailed survey of the cosmology, 
emphasizing the novel aspects that are a result of the D3 
symmetry. 

An immediate consequence of the non-Abelian symme- 
try is that since there are two stable DM candidates for 
Tvirf-j < 3m x , the individual relic abundances are not fixed 
by the WMAP measurement of the DM relic density, but 



rather the sum is fixed: 



fl v h 2 



+ ^x+x*h 



O WMAP 1,2 
"dm 11 



0.1126±0.0036 [28|. This means that different choices for 
a certain combination of parameters will fit the WMAP 
measurement while shifting the fractional abundances of 
77 and x- 

Annihilation into SM. The first class of reactions (a) 
in which 77 and x annihilate into light SM particles is 
common to any scalar DM model with a Higgs portal 
interaction. In the case m x ^ n <C m^, the thermally aver- 
aged annihilation cross section may be written as 15 - 17j | 



Aajv 



2„,2 



(Am 2 



m 



I) 2 + mlTl 



(10) 



where i = 1 (2) corresponds to 77 (x) and Xsm de- 
notes all kinematically allowed final states (e.g., Xsm = 
ZZ, WW, bb . . . for m, > mz)\ v in the nominator of the 
right hand side stands of course for the Higgs expectation 
value in Q. Furthermore, Th denotes the decay width 
of the SM Higgs boson (including decays to the new dark 
scalars if kinematically accessible) and we have defined 
r, = rV-vXsAf ( m h* = 2m<) to denote the decay width 
of a virtual SM Higgs of mass m^* — 2mj. For scalars 
above the weak scale, m,^ 3> m^, the annihilation cross 
section becomes 



Airmf 



(11) 



For the lightest of the DM species, annihilation into 
SM final states is particularly important since the DM 
conversion processes ( c ) are likely suppressed by kine- 
matics and Boltzmann statistics. Consider, for concrete- 
ness, the case m v < m x for which r\ dominantly depletes 
via annihilation into the SM final states. This implies 
a sizable Higgs portal coupling ol\ in order to obtain an 
acceptable relic density. For example, for <C nih and 
below the s— channel resonance at = m^/2, the dom- 
inant annihilation channel is r\r\ — > bb, with cross section 
(av) --y^g — Zu\m 2 / Tim^. One then finds a minimal 
coupling otx > 0.2 x (m/j/120 GeV) 2 . Also in the oppo- 
site case, when 77 is heavier than the SM Higgs h (but 
still lighter than x)i the annihilation cross section (fTT|) 
dictates a minimal coupling ot\ > 0.1 x (m n /600 GeV). 
Only near resonance, m n = mh/2, can the coupling ol\ be 
much smaller while still allowing for sufficient depletion 
of 77 particles. Conversely, the heavier species x need not 
couple strongly to the Higgs portal since the DM conver- 
sion processes (c) are efficient for m x > m^. 

Semi- annihilation. Consider now the case when x is 
lighter than 77. The Higgs portal coupling 0.2 should now 
be sizable to give x the right relic abundance. In addi- 
tion, the semi-annihilation processes (c) become available 
for m x > rrih . Indeed, this may be the dominant process 
for which one can obtain a semi-analytic solution to the 
relic abundance of x 22| : 



n x+x ,h 2 ~2x 1/2 



1.07 x 10 9 a; F 



(12) 



3* denotes the effective number of relativistic degrees of 
freedom at the time of the freeze-out, and the Planck 
scale is M P = 1.22 x 10 19 GeV. The decoupling tem- 
perature Tp can be determined upon solution of xp ~ 
log{0.038c(c + l)[±((Tv) xx ^ x . h ]m x Mpg- 1/2 x F 1/2 } with 
c = v2— 1, where x — m x /T. In the D3 model the ther- 
mally averaged semi-annihilation cross section is given 
by ' 



(av) 



ia\^i\v 2 
32ttto6 



XX^hx* 



h L h 



(13) 
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FIG. 1. Semi- annihilation: The comoving density Y(x) vs. 
x = m x /T for \- The thin dotted line shows the equilibrium 
value Y eq (x), while the thick lines show the evolution of the 
comoving density. We take m x = 200 GeV, = 120 GeV, 
ci2 = 0.016, and ot\ — Q3 = Q4 = while assuming m v S> m x . 
With y,\ = (dashed) \ freezes out too early and overdoses 
the universe whereas with fi\ = 1000 GeV (solid) \ freezes out 



later so that the relic abundance is Q 



x+x' 



such a small value of the Higgs portal coupling a.%. This 
illustrates that semi-annihilation allows for efficient depletion 
of x if the cubic coupling fii is sizable when m x > ra/,. 



We see from Eq. (|T5j) that both the Higgs portal cou- 
pling a2 and the cubic coupling fj,i must be nonzero for 
semi-annihilation to play a role. Because 6*2 ^ 0, anni- 
hilation into SM final states (a) can occur as well and 
will in general dominate unless the coupling fix is rela- 
tively large. Furthermore, because of the steep falloff of 
the semi-annihilation cross section (| 13[) with increasing 
rn x , the process is most effective when x is n °t much 
heavier than the SM Higgs, m x > rut- The effect of 
semi-annihilation is illustrated in Fig. [I] which shows the 
evolution of the comoving density Y = n/s, where s is 
the entropy density. 

DM Conversion. We next consider the processes of 
DM conversion (d), which consist of 2 — > 2 reactions 
that convert DM of one species into another. There are 
two types of reactions within class (d) that may be dis- 
tinguished. The reactions in the first line of (d) may be 
regarded as "dark" semi-annihilation (e.g., XX ~> VX*i 
etc.), in which the number of DM particles of a particu- 
lar species changes by one unit. Such processes are gov- 
erned by the nontrivial interaction predicted by the D3 
symmetry in Eq. ([6]) with coupling Q4. The reactions in 
the second line of (d) are regarded as "dark" annihilation 
(e.g., XX* ~^ 7 1 7 1 an d the inverse process), which change 



the number of each dark species by two units. The dark 
annihilation processes can be mediated by the Higgs por- 
tal when ai 2 are sizable, or via the direct scalar coupling 
a 3 in Eq. ©. 

The DM conversion processes that reduce the number 
of the heavier DM species, be it n or Xi are kmemati- 
cally favorable and hence very efficient at depleting such 
heavy particles. Thus DM conversion may have a dra- 
matic influence on the eventual number density of the 
heavier species. Unless the masses of DM particles are 
nearly degenerate, the relic abundance of the lighter DM 
species will not be influenced. 

Let us illustrate the effect of DM conversion by assum- 
ing for the moment that m x 3> m„ for which there are 
two stable DM species. As discussed above, 77 — being the 
lighter state — must have a sizable Higgs portal coupling 
oil. On the other hand, x may deplete via DM conversion 
processes. For example, if the dominant process is XX ~^ 
tix*} the relic abundance may be estimated from (|12[) 
with the replacement (o~v) xx -), x *h ~> ( av )xx^>vx* ■ ^ n our 
minimal D3 model the "dark" semi annihilation cross 
section is found to be 



3a 2 



((TV) 



I28irm 2 x ' 



(14) 



One needs in this case a minimal coupling «4 > 0.05 x 



despite (m x /50GeV). In contrast to semi-annihilation into the 



SM Higgs (flU)) . "dark" semi-annihilation ([T4"f is quite 
efficient and requires only a modest coupling on to sig- 
nificantly influence the relic abundance of x- We illus- 
trate the effect of DM conversion in Fig. ([2]), focusing on 
"dark" semi-annihilation. A similar discussion applies to 
"dark" annihilation (e.g., XX* ~^ Wi etc). 

Late decay of 77. The final novel aspect of the cos- 
mology in the minimal -D3 model is the possibility of the 
late decay of 77. If m v > 3m x , the state 77 will decay 
via 77 — ¥ 3%, 3x* due to the coupling 04. In the limit 
3> 3m x , the total width of 77 is approximately given 

l>v 



11 ~ 1536tt 3 
and which implies a lifetime of 

/10" 



10" 



V a 4 



(15) 



HM] . (16) 



For small values of (which are technically natural) , 77 
becomes long-lived in the sense that it decays after chem- 
ical decoupling, 3> ^,f- The cosmic time at freeze-out 
reads t n>F ^ lCT 8 s x (m r) /100 GeV)" 2 (^ :F /20) 2 , where 
x n = m v /T and T Vt p is the decoupling temperature. For 
a given , it is then easy to read off the maximum values 
of oa from (I16|) so that the late-decay condition t v t n .F 
holds. Clearly, the decays 77 — > 3%, 3%* can then sig- 
nificantly contribute to the relic abundance ^ x + x * h 2 — 
provided that ^-annihilation has frozen out. In the limit 
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FIG. 2. DM Conversion: The comoving density Y(x) vs. 
x — m x /T for n (gray) and \ (black). The thin dotted lines 
show the equilibrium value Y eq (x), while the thick lines show 
the evolution of the comoving density. We take m v = 5 
GeV, m x = 50 GeV, m h = 120 GeV, ai = 0.45, and 
Mi = Q.2 — «3 = 0. With on = 0.01 (dashed black) x 
freezes out too early and overdoses the universe whereas with 
a.4 = 0.075 (solid black) \ freezes out later so that the total 
relic abundance is Q v + fi x +x* ~ f2i3M AP ■ This demonstrates 
that DM conversion is very efficient at depleting the heavier 
DM species \- Conversely, the relic abundance of the lighter 
species 77 is the same for both values of 04 and therefore not 
affected by DM conversion. 

where \ is initially depleted, late 77-decays lead to a non- 
thermal ^-population with relic DM abundance 

fl x+x *h 2 =3^Q v h 2 , (17) 

where fl n h 2 represents the would-be relic abundance of 
77 particles today, had they not decayed. The late decay 
scenario is illustrated in Fig. [3] 

For t v > Is, the late-decay scenario (d) is cosmologi- 
cally constrained. The main decay mode 77 — > 3%W is in- 
evitably accompanied by 77 — > 3x^+h* where the virtual 
h* decays into kinematically accessible final states. This 
leads to electromagnetic and hadronic energy injection 
into the primordial plasma. If the decay happens dur- 
ing or after Big Bang nucleosynthesis (BBN) the induced 
spallation reactions of the light elements can spoil the 
concordance of BBN abundance predictions with their 
observationally inferred values (for a review see, e.g., 
[29| .) Since energy depositions as small as O(fewMeV) 
per baryon can probed in BBN [30], branching fractions 
of the 3% + h* decay mode of, say, O(10~ 4 ) are still 
constrained. To circumvent this issue as well as poten- 
tial warm DM constraints due to the free streaming of 



FIG. 3. Late decay: The comoving density Y(x) vs. x = 
m x /T for 77 (gray) and x (black). The thin dotted lines show 
the equilibrium value Y eq (x), while the thick lines show the 
evolution of the comoving density. We take m n = 800 GeV, 
m x = 200 GeV, m h = 120 GeV, ati = 0.11, a 2 = 0.15, 
c*4 = 10 -9 and /xi = 0:3 = 0, leading to a long lifetime t v ~ 
40/xs. Initially x is depleted through efficient annihilation into 
SM particles, whereas 77 freezes out with a significant density. 
Later, 77 decays and repopulates x> X* which then comprise 
the observed DM. 

the decay products x, one can simply impose r, ; < 1 s. 
For a kinematically unsuppressed decay, this requires 
|a 4 | > lO^V 100 GeV/mr,. 

Finally, we note that in the limit cu — > 0, 77 becomes 
stable even for m ri > 3m x . This is easily understood 
since in this limit the Lagrangian displays a Zi symmetry 
j] — > —77 which renders 77 stable. 

Direct detection phenomenology. The minimal 
D3 model of DM can be efficiently probed by direct detec- 
tion experiments. Furthermore, compared to the canoni- 
cal Z2 scalar singlet model, the non-Abelian D3 model of- 
fers several distinct direct detection signatures which can 
be observed by current and next generation experiments. 
We now examine the constraints and future prospects 
offered by these experiments. 

We use the standard formalism [3l| to predict the nu- 
clear recoil rates from elastic scattering of 77 and x ■ The 
differential event rate is given by 

£t=Nt*-[ *vvf(v,v.)£L, (18) 

where i = 1(2) corresponds to r}(x)i Nt is the number of 
target nuclei per unit detector mass, pi is the local DM 
mass density for a given species, and f(v,v ) is the DM 
velocity distribution taken to be Maxwellian and trun- 
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cated at an escape velocity of v esc = \v + v c \ = 600 km/s; 
v m in is the minimum velocity required to cause a nuclear 
recoil with energy Er and v e and v are the velocity of the 
earth in the galactic frame and the velocity of the DM 
particle in the earth's frame, respectively. The theoreti- 
cal rate in Eq. (fT5)l is corrected to account for potential 
quenching of the recoil signal as well as for finite detector 
resolution, efficiency, and acceptance to obtain a predic- 
tion for the observed rate. The total number of events in 
a certain energy interval is then found by integrating the 
rate and multiplying by the raw exposure. 

The microscopic physics of our D$ model enters into 
the scattering cross section in Eq. (TT51) . The scalar rep- 
resentation of Z?3 gives rise to spin-independent DM- 
nucleus scatterings. The differential cross section is con- 
ventionally expressed as 



dE h 



m N a n 



2v 2 



f p Z + f n (A-Z) 



fn 



F 2 (E R ), (19) 



where tojv is the mass of the target nucleus, p n is the 
DM-nucleon reduced mass, and F 2 (Er) is the nuclear 
form factor; we take the Helm form factor following [31] 
with nuclear skin thickness s = 0.9 fm. Spin-independent 
elastic scattering of the DM particles r\ and \ with nuclei 
is mediated via a t— channel exchange of the SM Higgs 
boson and therefore governed by the couplings ot\ and 
a2, respectively. The DM-nucleon cross section is given 
by 



(20) 



where /„ = £\ 



DM-nucleon coupling 32]. 



n) ~ 0.52 GeV is the effective 



Because there are two species of DM in our model, 
the local density of each species is typically not equal to 
po ~ 0.3 GeV cm -3 . As in Rcf. [33] we assume that the 
local density of each DM particle is proportional to its 
cosmological abundance, i.e., 



oWMAP 
1 l DM 



PO- 



(21) 



For a given set of model parameters, we compute the relic 
density fij of each species by integrating the Boltzmann 
equations, find the local DM mass density pi according to 
Eq. (|2"Tj) . and calculate the nuclear scattering rate (fT8|) . 

Experiments. In the following discussion we shall fo- 
cus on the liquid xenon experiments in anticipation of the 
upcoming XENON100 one-year data release as well as 
the prospect of future ton-scale experiments which have 
the potential to probe large regions of the D 3 model 
parameter space. Furthermore, past data sets from 
XENON10 [al El] and XENON100 [H] yield current 



limits on spin-independent scattering that are competi- 
tive with other experiments, such as CDMS II (37l. l38j|. 



The XENON experiments use a prompt scintillation light 
signal (SI) and a delayed ionization signal (S2) to dis- 
criminate between nuclear and electron recoil events. Re- 
cently, there has been a lot of debate on the scintil- 
lation yield for nuclear recoils, C c g, as well as on the 
stochasticity of the SI signal for lowest Er. Since £ e ff 
defines the detector threshold, this quantity is of im- 
portance when inferring cr„-limits for small DM masses 
O(10 GeV) [H EJ]. We use the results of the detailed 
study [13] to account for resolution, efficiency, and accep- 
tance of the XENON10 H and XENON100 detec- 
tors where we employ the most conservative assumption 
on £ G ff (yielding the weakest limits on o~ n .) Constraints 
are obtained by Yellin's maximum gap method 41] which 
accounts for the O(10) events observed in the XENON10 
re-analysis [35| of its 316kgxdays data sample. The 
XENON100 collaboration has published an analysis of an 
11 day run yielding a raw exposure of 447 kg x days [36j |: 
no candidate events were observed. 

In addition, a data release from a one year run of 
XENON100 is expected shortly. To estimate the reach of 
this data set we assume an exposure of 100 live days, col- 
lected with their fiducial detector mass of 30 kg. Lastly, 
for the sake of exploring future sensitivit y of a ton-scale 
liquid xenon experiment, following Ref. [42j j we assume 
less than one background event; we take a (moderate) 
raw exposure of ltonxyr. We mimic energy resolution 
and efficiency of the detector by employing the same 
analysis [40( as used for the XENON10 detector above; 
this implies a more optimistic assumption on the low re- 
coil acceptance then had we used the similar analysis for 
XENON100. 

As we now show, large portions of the parameter space 
in the D3 model have been or will be probed by liquid 
xenon experiments. We will describe the constraints and 
future sensitivities in detail for two novel scenarios that 
may occur in the D3 model: 1) the late-decay scenario 
(m n > 3to x ), and 2) the two-component DM scenario 
(m n < 3m x ). 

Constraints for m v > 3m x . If 77 is more than three 
times as heavy as x, it is unstable, dictating that \ ac ~ 
counts for all of the cosmological DM. Without making 
assumptions about the lifetime of 77 (or equivalently the 
coupling 0^4), the Higgs portal coupling a-i can take a 
range of values and still be consistent with the relic abun- 
dance constraint. The minimum allowed value of a2 is 
determined by the requirement that x does not overdose 
the universe upon freezeout. However ct2 can in fact be 
much larger than this minimum value because late de- 
cays of rj can replenish an initially depleted population 
of x particles. 

In Fig. [4] we show the constraints in the m x -ai plane 
for the case m v > 3m x and p = (so that semi- 
annihilation XX X*h is absent). Along the thick solid 
line the relic abundance of x is m concordance with that 
inferred from observation, S! x + X * = S!^m Ap . The re ~ 
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FIG. 4. 

plane for m v > 3m x and /i = 0. Along the thick solid (black) 
line fi x +x* = fi™ AP holds. The region below is excluded, 
f2 x + x * > f^^^J Ap . Only x ls stable so that above the solid 
black line n — ► 3x^ populates x'*' such that p x = po- The 
gray shaded region is excluded by current direct detection 
experiments. As labeled, the lines depict the individual ex- 
periments and show their current (projected) sensitivity for 
exclusion. For m x ~ rrih/2 = 60 GeV primordial annihilation 
can proceed resonantly via Higgs exchange. The horizontal 
dotted line shows a vacuum stability constraint for A3 = 2; 
values above are excluded. 



gion below is excluded since ^ x + x » > ^dm ■ ^ n tne 
region above, the decay rj must source the relic 

abundance to obtain Sl x + X * = ^5m AP and thus a vi- 
able cosmology. We see that, except near the resonance, 
2m x ~ mh = 120 GeV, the Higgs portal coupling ct2 
must be sizable. The shaded region is excluded due to 
measurements of past direct detection experiments. In 
particular, the light DM region with m x < 10 GeV is 
most constrained by XENON10 (as well as by the CDMS- 
II low threshold analysis [H| ) , while heavier DM is most 
constrained by the XENON100 first run with 11 live 
days (as well as the final exposure of CDMS II [37|). 
We also show the projections for the upcoming release 
of XENON100 one-year data, as well as a future ton- 
scale xenon experiment. While the former is not ex- 
pected to improve upon the XENON10 low-mass limit, 
the latter will not only improve upon this limit but will 
essentially probe all values of m x that lie to the right 
of the Higgs resonance region. Finally, the dotted hor- 
izontal line shows the requirement that the electroweak 
vacuum JJJ) is a global minimum; values of a.i above that 
line will make the vacuum configuration metastable. The 
bound depends on the quartic coupling A3 (here A3 = 2) 
and also on the sign of ui. 



FIG. 5. Direct detection constraints: Shown is the m v -ai 
plane for m v < m x . Along (below) the thick solid line Q v = 
Qom AP (M n > ^dm AP ) holds. The thin gray lines parallel 
and above are contours with a fraction of 0.1 and 0.01 of 
C7 (bottom to top.) The near-to vertical lines depict 
the individual experiments and show their current (projected) 
sensitivity for exclusion. The Higgs resonance region is shown 
more closely in the inlay. 



Constraints for m v < 3m x . If m. n < 3m x , there are 
two DM species since the decay n — > is kinemati- 

cally forbidden. The individual relic abundances of rj and 
are not constrained by WMAP, and so it may appear 
that further model assumptions are necessary to con- 
strain the two-component scenario. In fact, robust direct 
detection constraints exist for the lightest DM species 
independent of the fractional composition of the cosmo- 
logical dark matter. For the sake of discussion, let us 
assume that 77 is the lightest DM particle, m v < m x . 
As discussed in detail above, the only kinematically fa- 
vorable annihilation channel for 77 is into the light SM 
states (jlip . This means that r\ must have a minimum 
Higgs portal coupling ol\. Because the same coupling ot\ 
controls the DM-nucleon scattering cross section ([20)) , di- 
rect detection experiments constrain large portions of the 
m^-ai parameter space. 

In Fig. [5] we show constraints in the m^-ai plane for 
the case m v < m x . The lines as labeled are analogous 
to the ones shown in the previous Fig. |H This time, 
however, the regions between 6 GeV < m v < 50 GeV 



and 70 GeV < m n < 80 GeV are excluded for all val- 
ues of a\. This may seem counterintuitive at first sight, 
since as o.\ increases the relic abundance decreases: f2, ; ~ 
1/ (pv) m -+x S M x a \ 2 so that Pr) cx aj~ 2 by Eq. I|21[) . 
This can also be seen by the contour lines parallel to the 
thick solid line which now show decreasing fractions 0.1 
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and 0.01 (from bottom to top) of the DM abundance. 
For increasing ct\ there is, however, a corresponding in- 
crease in the DM-nucleon scattering cross section (|2U|) : 
all x a x- Since dR^/dEn cx p^cr^ the dependence on 
at\ cancels out in the nuclear scattering rate (fT51) . This 
explains why the various direct detection constraints as 
labeled are now near-to vertical lines. 

While the constraints on the lightest DM particle (as- 
sumed to be r\ for the discussion) are not strongly depen- 
dent on the assumptions made about other model pa- 
rameters, the same is not true for the heavier species 
(here assumed to be \). This is because the coupling a 2 
governing the interactions of \ with the SM is largely un- 
constrained by the cosmology, since x may efficiently de- 
plete through DM conversion processes ( e.g., XX ~^ X* 7 ?; 
etc.). Hence further assumptions must be made about 
other couplings in the Lagrangian (|6]) to make concrete 
statements about the heavier DM component. 

Discovering two-component dark matter. The two- 
component nature of DM for m v < 3m x is one of 
the most intriguing features of the non-Abelian D3 dis- 
crete symmetry. Naturally, the question arises: can 
one discriminate between a single-component and a two- 
component DM scenario at future direct detection exper- 
iments? Here we will attempt to provide a quantitative 
answer to this question. Our analysis will be based on 
the shape of the nuclear recoil spectrum, which may ex- 
hibit distinct features depending on whether one or two 
components are scattering in the detector. 

To answer the question posed above we Monte Carlo 
generate a sample recoil spectrum for a ton-scale liquid 
xenon experiment with m v = 5 GeV and m x — 200 GcV. 
We further set ct\ — 0.45, ct2 — 0.065 and 0:4 = 0.3 for 
which the dark matter is distributed equally, fl v = Q x = 
^Sm AP /2 (note that DM conversion is used to deplete 
x). This parameter point is consistent with current di- 
rect detection constraints (see Fig. [5]). Since the XENON 
experiments have a poor energy resolution, the bin-width 
is chosen to be 4keV r with a total number of ten bins, 
starting from 2keV r . The total number of events to be 
generated is then drawn from a Poisson distribution with 
its mean given by the number of theoretically expected 
events. A sample "experimental" spectrum is shown by 
the dots in Fig. [6l The events in the first bin are dom- 
inated by scatterings of the light state r/ on the Xe tar- 
get. In the second bin both states r\ and X contribute 
equally. At higher recoil energies all events are essentially 
due to X ■ This behavior can of course be traced back 
to the favorable relation <C m x which yields rather 
distinct exponential forms of dR/dE^. 

We now attempt to discriminate the two-component 
origin of the generated spectrum from a single- 
component DM hypothesis. For this we minimize the 
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FIG. 6. Recoil spectrum: The dots show a Monte Carlo 
generated recoil spectrum of the two-component DM scenario 
with input values m v = 5 GeV, ai = 0.45, m x = 200 GeV, 
c*2 = 0.065, and pi t % — po/2 for a ton-scale liquid Xe detector 
with raw exposure of 1 ton x yr. The dashed line (blue) shows 
the attempt to explain the data with a generic single DM 
particle. The solid line (black) is obtained from a DM model 
with two-components; see main text for details. 



Poisson log-likelihood function 



Xl = 2 



N. 



N ohs + N obs m (jv° bs /j\^ h )] (22) 



bins i 



to obtain the best fit and also to assess the goodness-of- 
the-fit Here, iV° bs (N% h ) is the number of observed 
(theoretically predicted) events in the i-th bin. We use 
the program MINUIT [3] to minimize ([22]). Note that 
with the procedure outlined above, we only account for 
statistical uncertainties in the event rate; a more detailed 
analysis incorporating systematic uncertainties goes well 
beyond the present scope. 

We first attempt to fit the generated recoil spectrum 
with a single DM particle. The DM mass and spin inde- 
pendent DM-nucleon cross section are used as the fitting 
parameters. Using the procedure described above, one 
obtains the dotted (blue) histogram in Fig. [6] As can 
already be seen by eye, the one-component fit is poor. 
Indeed, from the goodness-of-fit test the single-DM hy- 
pothesis is rejected with confidence well above 99%. 

We now turn to a two-component DM model. In this 
case the recoil spectrum depends on the DM masses nii, 
cross sections, o~n\ and local densities pi, with z = 1,2 
and subject to the constraint p\ + pi = pq. Thus, there 
are in principle five parameters to be inferred from the 
data of Fig. O There exist, however, degeneracies be- 



9 



tween different parameters in the two-component spec- 
trum so that only certain combinations are accessible. 
To understand this it is helpful to examine the paramet- 
ric dependence of the recoil spectrum (fTS)) on n%i, pi, 
and cin\ Using (|2"0|) . and in the limit of \v c \ — and 



dRi 



i one finds 

(0 
PiCTn 
ex. exp 

mi 



m N 
mi 



(23) 



where vq is the DM velocity dispersion. Observe that 
1) only the product pi<7„ enters in dR/dE R and 2) the 
exponential shape of the spectrum becomes independent 
of the DM mass for rrii S> toat. From the first observation 
it is clear that one cannot independently determine both 
Pi and On j but only the product. The second observation 

(i) 

suggests that a second degeneracy between piOn and m, 
develops once m; >mj(. Therefore, further assumptions 
about the model parameters are required if one desires to 
break these degeneracies. Fixing, for example, mass and 
density of the heavier DM to its input values, one then 
easily recovers recovers all remaining model parameters 
within a "ler" -range or better. The result of this fit is 
shown by the solid (black) histogram in Fig. |6l 

Note that the successful determination of the param- 
eters of the lighter state required that at least two bins 
were populated. If one tried to fit the light DM particle 
to a single bin, many combinations (mi, <J\) would repro- 
duce the signal thus leading to a high degree of degener- 
acy in mi and a±. In this regard, it is important to note 
that — complementary to liquid noble gas experiments — 
potential future ton-scale cr yog enic detectors such as su- 
perCDMS or EURECA [46| may be very powerful in 
disentangling the parameters of a multi-component DM 
scenario. The reason lies in a much better energy reso- 
lution. For example, Ref. [42[ estimates for a ton-scale 
Ge detector a Gc < 0.5 keV (E R < 40keV). This will al- 
low one to obtain a more detailed picture on the spectral 
shape of the signal. Moreover, a different target mass 
with respect to Xe may prove most valuable when at- 
tempting to draw differential conclusions by combining 
results from both detector designs. 

Discussion. In this paper we have investigated the 
simplest model in which DM is stabilized by a non- 
Abelian discrete symmetry. The model is based on the 
symmetry group D3, which is the smallest non-Abelian 
finite group. The non-Abelian nature of this theory man- 
ifests itself in the matter content and interactions of the 
dark sector, which leads to multi-component DM and a 
novel early universe cosmology. Robust constraints from 
direct detection experiments exist for the lightest DM 
species, and the two-component nature of DM can po- 
tentially be tested at future ton-scale experiments. 

While we have focused on the phenomenology relevant 
for direct detection experiments, the D3 model also has 
implications for high energy colliders and indirect de- 



tection probes. At colliders, the most important con- 
sequence of the D3 model, like other models in which 
DM couples to the SM through the Higgs portal, is to 
cause the Higgs boson to decay invisibly if kinematically 
allowed. An invisibly decaying Higgs boson can be dis- 
covered using the vector boson fusion and ZH produc- 
tion channels at the LHC (y/s = 14 TeV) with 30 ftr 1 
of data [47j | . Another generic DM collider probe is a 
monojet produced in association with pairs of dark mat- 
ter particles 48, 49|. In scalar singlet DM models, this 
occurs via gluon fusion accompanied by emission of a 
jet. While the Tevatron is not sensitive to the monojet 
signature in this class of DM models, the LHC could po- 
tentially probe regions with large Higgs portal couplings 
a i,2 ~ 0(1) [i^]. We also note that a more complete 
study of the monojet signature in Higgs portal DM mod- 
els going beyond the effective field theory approach is 
warranted given that a relatively light Higgs and a top 
quark loop enter into the amplitude. 

A potential indirect probe of the D3 model is to ob- 
serve gamma rays from annihilating rj or % particles. In 
particular, Fermi-LAT data on the isotropic gamma-ray 
diffuse emission [5(|, depending on various astrophysi- 
cal assumptions, can constrain low mass < 10 GeV DM 
particles annihilating through the Higgs portal [5l|. For 
such constraints to be relevant, the lightest component 
should carry most of the cosmic DM abundance, since 
the gamma-ray flux scales as the square of the DM den- 
sity. In particular, light DM particles with very large 
Higgs portal couplings will be depleted efficiently and 
therefore not subject to such constraints. Another possi- 
bility is to observe monochromatic gamma-ray lines via 
VVi XX* ~ ► h* — > 77. This has been studied in Ref. |52j |. 
and Fermi-LAT data allows one to probe the model deep 
within the resonance region m, hx ~ m^/2. 

There are many directions for future investigations 
with DM and non-Abelian discrete symmetries. Within 
the minimal D3 model it would be interesting to con- 
sider the case in which r\ (x) condense, leading to the 
spontaneous breaking of the -D3 symmetry to the Z3 (Z2) 
subgroup. This will still lead to a viable model of DM 
and likely will have distinct signatures compared to the 
unbroken D3 symmetry. It would also be interesting to 
explore larger non-Abelian finite groups, and ultimately 
non-Abelian discrete gauge symmetries as the origin of 
DM stability. Finally, given that non-Abelian discrete 
symmetries may underlie the patterns observed in the 
quark and lepton masses and mixings, it is interesting to 
speculate that DM stability may intimately be connected 
to flavor symmetries. 

Stability on cosmological time scales constitutes one of 
the few and robust guiding principles in formulation of 
a theory of dark matter. Alternatives to the canonical 
Z2 parity symmetry may lead to new phenomena associ- 
ated with dark matter, as has been clearly demonstrated 
by our study of the minimal model of D3 dark matter. 
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We anticipate that new data from a variety of experi- 
mental and observational fronts will soon bring us closer 
in unraveling the particle nature of dark matter and the 
symmetry which protects its lifetime. 
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A Boltzmann Equations 

The Boltzmann equation describing the evolution of 
the number density is given by 

dn ■ 

-r£+3Hn i = C i) (A.l) 

where i = rj, x and H(T) is the Hubble parameter ac- 
counting for dilution of rii due to the cosmological ex- 
pansion of the Universe. The collision terms Ci have the 
following contributions: 
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xx" ^nn 
(A.3) 



If the reheating temperature of the radiation domi- 
nated Universe was high enough, all particles in con- 
sideration were once in thermal equilibrium. Assuming 
Maxwcll-Boltzmann statistics and instant thermalization 
of the reaction products, Ci can be written in familiar, 
integrated form. Furthermore, CP invariance relates var- 
ious collision terms so that we find for the individual con- 
tributions of C„: 



a 



a 



VV^XX 



a, 



»)X^X*X* 



iwn^xx 



IVX^X'X* 



2 « q ) 

2 



cq '"n 



C n 



= c 



vx->x x ' 



Cn 



= -(T 



V^3X/ 



("x q ) 3 "" q 



c 



»?-S-X*X*X* 



C n 



(A.4) 



while the terms in C x in Eq. (IA.3I) are written as 
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(A.5) 



The thermally averaged annihilation cross section for 
process a + 6 — > c + d is in general given by 
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where the quantity W ab is defined as 
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In W ab the squared matrix element \M. ab ^cd\ 2 is inte- 
grated over the cm. scattering angle 9. Ki denotes the 
modified Bessel function of order i and p ab is the mo- 
mentum of relative motion of a and b. Furthermore, g a , 
g b are the number of internal degrees of freedom, S cd is 
a symmetry factor accounting for identical final states, 
and X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. For the 
decay a — > b + c + . . . , the thermal average reads 



/r v _ Ki(m a /T) r 

K 2 (m a /T) 



(A., 



where r a ^fc+ c +... is the standard decay rate for a particle 
a at rest. 
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